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We review the detectability of gravitational waves generated by oscillations excited during a 
phase transition from hadronic matter to deconfined quark-gluon matter in the core of a neutron 
star. Neutron star properties were computed using a Boguta and Bodmer's based model and the 
MIT bag model. The maximum energy available to excite mechanical oscillations into the star is 
estimated by energy difference between the configurations with and without a quark-gluon matter 
core. On basis of the planned sensitivity of present laser interferometers (VIRGO or LIGO I) and 
those of the next generation (LIGO II), the maximum volume to be probed by these experiments is 
determined. These results are used as an indication of the potential detectability of neutron stars 
as sources of gravitational waves. Our results indicate that the maximum distance probed by the 
detectors of the first generation is well beyond M31, whereas the second generation detectors will 
probably see phase transition events at distances two times longer, but certainly not yet attaining 
the Virgo cluster. 
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I. INTRODUCTION 

The first generation of large gravitational interferomet- 
ric detectors as the French-Italian VIRGO and the Amer- 
ican LIGO, should be fully operational within one or two 
years. The best signal-to-noise (S/N) ratio that can be 
achieved from these detectors implies the use of matched- 
filter techniques, that require a priori the knowledge of 
the signal waveform. Thus, the identification of possible 
sources having a well defined signal is a relevant problem 
in the detection strategy. 

Neutron stars are certainly one of the most popular po- 
tential sources of gravitational waves (GWs), since they 
can emit by different mechanisms having a well known 
waveform. Rotating neutron stars may have a time- 
varying quadrupole moment and hence radiate GWs, by 
either having a triaxial shape or a misalignment between 
the symmetry and the spin axes, which produces a pre- 
cessional motion^ ||, 9. Moreover, fast rotating proto- 
neutron stars may develop different instabilities such 
as the so-called Chandrasekhar-Friedman-Schutz (CFS) 
instability^ ||, responsible for the excitation of density 
waves traveling around the star in the sense opposite to 
its rotation, or undergo a transition from axi-symmetric 
to triaxial shapes through the dynamical "bar-mode" 
instability^. All these mechanisms are potentially able 
to generate large amounts of energy in the form of GWs. 

Nonradial oscillations are also a possible mechanism 
for neutron stars emit GWs and this possibility was al- 
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ready discussed in the late sixties |7)]. One of the diffi- 
culties with this mechanism concerns the energy source 
necessary to excite the oscillations. The elastic energy 
stored in the crust and released by tectonic activity was 
recently considered ||, but the maximum available energy 
is likely to be of the order of lfj 44-45 erg. Thus, even if 
all this energy could be converted into nonradial modes, 
the maximum distance that a signal could be seen by a 
laser interferometer like VIRGO is only about 3.0 kpc[|). 
This result depends, of course, on the adopted equation 
of state, which fixes the mode frequencies and damping 
timescales. A considerable amount of energy would be 
available if the neutron star undergo a phase transition 
in the core. For instance, if quark deconfinement occurs, 
then the star will suffer a micro-collapse since the equa- 
tion of state of the quark matter is softer than that of 
the hadronic matter and the new equilibrium configura- 
tion will be more compact, having a larger binding en- 
ergy. The energy difference will partially cover the cost of 
the phase transition and will partially be used to excite 
mechanical modes, whose energy will be dissipated by 
different channels. The physical conditions of matter at 
densities above the saturation density (po = 0.15/m~ 3 ) 
is badly known and many first or second order phase tran- 
sitions have been speculated upon. The possibility of a 
quark and hadron mixed phase was considered in rcf . Q , 
where the conservation of the electric and the baryonic 
charges are satisfied not locally but within a volume in- 
cluding many droplets of quarks and of hadrons. This 
mixed phase does not include local surface and Coulomb 
energies of quarks and nuclear matter and, as a conse- 
quence of these terms, if the interface tension between 
quark and nuclear matter is too large, the mixed phase 
is not favored energetically jl0| |ll[ . In this case, the neu- 
tron star will then have a core of pure quark matter and 
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a mantle of hadrons surrounding it. In the absence of 
dynamical models able to follow the evolution of the star 
during the phase transition, we assume as in ref.[[l2|, that 
the structural rearrangement suffered by a "cold" star 
occurs on a dynamical timescale, i.e., of the order of mil- 
liseconds, shorter than the smooth global stellar readjust- 
ment suggested in refs.|l3|, [l4| and the gradual transition 
expected to occur in hot proto-neutron stars j^. Radial 
modes excited by the micro-collapse do not radiate GWs, 
but an important coupling with rotation exists Jl6[ and, 
if the star rotates, the oscillations will be damped not 
only by dissipation of the mechanical energy into heat 
but also by the emission of GWs. 

In the present work, the detectability of GWs gen- 
erated by oscillations excited during a phase transition 
in the core of a neutron star is reviewed p7[. Neutron 
star properties were computed using a description of 
the hadronic matter based on the work of Boguta and 
Bodmerp^], including the fundamental baryon octet, the 
isovector meson g and lepton degrees of freedom (see 
references [15] for details). Hybrid models, including 
a quark-gluon core, with the same baryonic number as 
the pure hadron configuration were also computed and 
the energy difference between both states was used as an 
estimate of the maximum energy available to either ex- 
cite mechanical oscillations or to be converted into heat. 
Then, using the planned sensitivity of present laser inter- 
ferometers like VIRGO (or LIGO I) and those of the next 
generation (LIGO II), the maximum volume of space that 
can be probed by these experiments is calculated, since 
this indicates the potential detectability of these sources. 
The plan of this paper is the following: in section 2 the 
star models are reviewed, in section 3 the gravitational 
wave emission is discussed and finally, in section 4, con- 
clusions are given. 
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This equation represents nuclear matter as composed by 
a mixture of the fundamental baryon octet (p, n, A, E + , 
E°, £~, E - , S°) coupled to three mesons (a,ui,g) and 
lcptons (for the details see [15]). The scalar and vec- 
tor coupling constants, g CT , g w and the coefficients b, c 
were determined by imposing that the model bulk prop- 
erties should be able to reproduce the binding energy 
Eb (= -16.3 MeV), the compression modulus K (= 240 
MeV) and the nucleon effective mass M* = M — g a a (= 
732 MeV) at the saturation density p (=0.153 far 3 ). 
Additionally, the isovector coupling constant g g is de- 
termined from the coefficient for the symmetry energy 
in nuclear matter, &4 (= 32.5 MeV). We have used the 
universal hyperon- nucleon coupling ratios Xi = 9 Hi 1 9%-, 
with i = cr,g,uj. The resulting coefficients used in our 
computations, are given in table 1. Figure 1 shows the 
equation of state derived from our model and in figure 2 
it is shown the energy density profile inside the star for 
two configurations having gravitational masses equal to 
1.2 and 1.6 M Q respectively. 



II. THE NEUTRON STAR MODEL 



The present model for the nuclear matter was already 
discussed by Marranghello et al.Q, who have investi- 
gated the effects of a finite temperature in the equation 
of state. Here, the same description for dense matter is 
adopted but only models with T=0 will be considered. 
Moreover, the coupling constants are slightly modified 
with respect to those considered in that work, in order 
to allow the maximum mass of the configuration to be 
compatible with recent results on the binary X-ray sys- 
tem Vela X-l[|o). For the sake of completeness, we recall 
here the main points of the model. 

The lagrangian density describing the nuclear matter 

is 



FIG. 1: Equation of state for hadronic matter (solid line) and 
for the quark-gluon matter (dotted line). 

Hybrid models including a quark-gluon core were also 
computed. The quark matter was described by the MIT 
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many parameters to be adjusted, even in the quark phase 
and specially in the hadron phase, we believe that only 
qualitative results can be obtained with some insights on 
the quantitative results. 



TABLE II: Physical conditions at the phase transition: energy 
densities and pressure are given in GeV fm -3 
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FIG. 2: Neutron star energy density-radius relation for 
M=1.6M© (solid line) and M=1.2M (dotted line) on which 
the quark-gluon core extends up to 8.6 km and 7 km, respec- 
tively. 



FIG. 3: Neutron star mass as a function of central energy 
density for hybrid star with constant pressure transition (solid 
line), hybrid (dotted line) and pure hadronis star (dashed 
line), configurations 



Lagrangian|2l| and the physical conditions at the decon- 
finement transition were estimated from the Gibbs cri- 
teria, namely, by the equality of the chemical potential 
and pressure of both phases, under conservation of the 
baryon number and electrical charge. The physical pa- 
rameters at the transition point are given in table 2, for a 
bag constant B 1 ^ 4 = !80MeV and a strange quark mass 
m s = 150 MeV. The gravitational mass of the star as a 
function of the central energy density is shown in figure 
3 for the case of a pure hadronic configuration (dashed 
line), a hybrid star with a quark core (solid line) and, 
for comparison, a case where a mixed phase exist in the 
center, compute in the same way as ref. H. 

It should be emphasized that our equation of state is 
quite steep ( ~ 2.6) near saturation and, as a conse- 
quence, the deconfinement transition occurs at densities 
just above the saturation value (p ~ 1.7po), producing 
hybrid stars with very extended quark-gluon cores. Here 
we take the opportunity to reiterate that, from the ac- 
tual status of theoretical predictions to the EOS, with so 



The maximum stable mass of a pure hadronic configu- 
ration is about M max » 2.1 M Q while for hybrid stars 
this limit is reduced to 1.73 M Q . Thus, the present 
calculations exclude the possibility that Vela X-l has 
a quark core. Notice that our models obey the Seidov 
criterium^U, namely, that the hybrid star will be sta- 
ble only if the energy jump across the transition surface 
satisfies the condition 



(2) 



where P, e q , e# are respectively the pressure, the energy 
density of quarks and hadrons at the transition point. 

One of the goals of this work is the determination of 
an upper limit for the energy able to excite the different 
oscillation modes of the star, when a modification in its 
internal structure occurs. In this sense, the simplest ap- 
proach is to compute the energy difference between two 
configurations having the same baryonic number, the first 
constituted of pure hadrons (H) , the second having a core 
of deconfined matter (HQ), and use such a difference as 
an indication of the maximum available energy. Table 3 
gives the parameters for five models defined by a given 
baryonic number |23j] 
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or the baryonic mass of the star (Mb ar = NMb); the 
expected gravitational mass of both configurations (pure 
hadron and hybrid), 
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expected radii, masses of the quark-gluon core and the 
energy difference obtained from the passage of configu- 
ration H to HQ, the binding energy. 



Eg = Mg - M bar . 



(5) 



III. THE GRAVITATIONAL WAVE EMISSION 

The transition from configuration H to HQ may occur 
through the formation of a mestastable core, built up by 
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TABLE III: Parameters of the stellar models: masses are in 
solar units and radii in km; H and HQ mean pure hadron and 
hybrid configurations respectively 
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M g (H) M g (HQ) M 3 (core) R (H) R(HQ) log (E) erg 



1.0749 
1.3202 
1.5724 
1.8342 
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1.0953 
1.2468 
1.4164 
1.6141 



0.1150 
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0.6013 
0.8626 
1.1766 



12.82 
13.04 
13.13 
13.01 
12.69 



12.47 
12.16 
11.98 
11.82 
11.50 



52.83 
53.28 
53.44 
53.52 
53.52 



an increasing central density. The increase in the cen- 
tral density may be a consequence of a continuous spin- 
down or other different mechanisms the star could suffer. 
This transition releases energy, exciting mainly the radial 
modes of the starjl7|. These modes do not emit GWs, 
unless when coupled with rotation jO^, a situation which 
will be assumed here. 

In order to simplify our analysis, we will consider that 
most of the mechanical energy is in the fundamental 
model. In this case, the gravitational strain amplitude 
can be written as 



h(t) = h e- {t/T ^- tulot) 



(6) 



where ho is the initial amplitude, ujq is the angular fre- 
quency of the mode and r gw is the corresponding damp- 
ing timescale. The initial amplitude is related to the 
total energy E g dissipated under the form of GWs by the 
relation 
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where G is the gravitational constant, c is the velocity of 
light and r is the distance to the source. 

Relativistic calculations of radial oscillations of neu- 
tron star with a quark core were recently performed by 
Sahu ct al.j2j|. However, the relativistic models com- 
puted by those authors do not have a surface of discon- 
tinuity where an energy jump occurs. Instead a mix- 
ing region was considered, where the charges (electric 
and baryonic) are conserved globally but not locally 
Oscillations of star models including an abrupt transi- 
tion between the mantle and the core were considered 
by Haenselp6| and Miniutti|27|. However a Newtonian 
treatment was adopted and the equation of state used in 
the calculations does not correspond to any specific nu- 
clear interaction model. In spite of these simplifications, 
these hybrid models suggest that rapid phase transitions, 
as that resulting from the formation of a pion condensate, 
proceed at the rate of strong interactions and affect sub- 
stantially the mode frequencies. However, the situation 
is quite different for slow phase transitions (the present 
case), where the mode frequencies are quite similar to 

In this case, scaling the 



those of a one-phase star 1 26 



results of p6| , the frequency of the fundamental mode 
(uncorrected for gravitational redshift) is given approxi- 
mately by 

(8) 



uq k, 63. 
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where the mass is given in solar units and the radius in 
km. 

Once the transition to quark-gluon matter occurs, the 
weak interaction processes for the quarks u, d and s 
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will take place. Since these reactions are relatively slow, 
they are not balanced while the oscillations last and 
thus, they dissipate mechanical energy into heat|2^]. Ac- 
cording to calculations of reference [21], the dissipation 
timescale can be estimated by the relation 
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where m s is the mass of the s-quark in MeV and M c is 
the mass of the deconfined core in solar masses. This 
equation is valid for temperatures in the range 10 8 — 
10 9 .K\ On the other hand, according to reference |ll| 
the damping timescale by GW emission is 
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where again the stellar mass is in solar units, the radius 
is in km and the rotation period P is in milliseconds. 

In a first approximation, the fraction of the mechanical 
energy which will be dissipated under the form of GWs 
is f g = ^ 1+t 1 . Notice that the damping timescale by 
GW emission depends strongly on the rotation period. 
Therefore one should expect that slow rotators will dissi- 
pate mostly of the mechanical energy into heat. In table 
4 is given for each star model the expected frequency of 
the fundamental mode (corrected for the gravitational 
redshift) , the critical rotation period (in ms) for having 
f 9 = 0.50, the GW damping for this critical period and 
the quality factor of the oscillation, Q =irvoT gw . 

After filtering the signal, the expected signal-to-noise 
ratio is 
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where h(v) is the Fourier transform of the signal and 
S n (u) is the noise power spectrum of the detector. Per- 
forming the required calculations, the S/N ratio can be 
written as 
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igw 



Q 2 



Sn{Vg W )J 1 + 4Q 



(14) 



5 



TABLE IV: Oscillation parameters: the damping timescale 
r gw is given for the critical period; maximum distances for 
VIRGO (V) and LIGO II (L) are in Mpc 
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In the equation above, the angle average on the beam 
factors of the detector were already performed |29|. 

From eqs.( 4) and (11), once the energy and the S/N 
ratio are fixed, one can estimate the maximum distance 
D mol to the source probed by the detector. In the last 
two columns of table 4 are given distances D maa; derived 
for a signal-to-noise ratio S/N = 2.0 and the sensitiv- 
ity curve of the laser beam interferometers VIRGO and 
LIGO II. In both cases, it was assumed that neutron stars 
underwent the transition having a rotation period equal 
to the critical value. 

We emphasize again that our calculations are based on 
the assumption that the deconfinement transition occurs 
in a dynamical timescale . In the scenario developed 
in ref . | fl3| , a mixed quark- hadron phase appears and the 
complete deconfinement of the core occurs according to a 
sequence of quasi-equilibrium states. The star contracts 
slowly, decreasing its inertia moment and increasing the 
its angular velocity until the final state be reached in 
a timescale of the order of 10 5 yr |L3|. Clearly, in this 
scenario no gravitational waves will be emitted and this 
could be a possibility to discriminate both evolutionary 
paths. 

IV. CONCLUSIONS 

The structure of pure hadronic configurations and that 
of hybrid stars having a quark core, with the same bary- 



onic number, were computed using a new equation of 
state The maximum stable mass for pure hadronic 
configurations satisfies the requirement of being higher 
than the mass of Vela X-l ( M = 1.86±0.16 M Q ), deter- 
mined recently |2(J from a new study of its orbital motion. 

As one should expect, the masses of quark cores in- 
crease with the baryon number of the configuration, as 
well as the difference between the binding energies be- 
tween hybrid and single phase objects. On the basis of 
this result and from the point of view of the GW emis- 
sion, it would be natural to expect that phase transitions 
occurring in more massive stars would be more easy to 
detect. In fact, this is not the case because, in the one 
hand the damping timescale due to the emission of GWs 
is inversely proportional to the mass and hence massive 
stars have low oscillation quality factors and, on the other 
hand the detectors have sensitivities optimized at fre- 
quencies around 200 Hz. In comparison with non-radial 
oscillations (for instance m= 1=2 modes), these have 
comparable frequencies but damping timescales one or- 
der of magnitude higher than those resulting from radial 
modes coupled to rotation. In this case, higher quality 
factors Q may be obtained but most of the mechanical 
energy will be dissipated into heat. In order to be an ef- 
ficient source of GWs by the mechanism here considered, 
the star must be a fast rotator, i.e., to have rotation peri- 
ods of the order of few milliseconds. Even in the most fa- 
vorable case (model 2, corresponding to a baryonic mass 
equal to 1.32 M©), if the rotation period is 4.0 ms, only 
1% of the available energy will be emitted as GWs. 

Inspection of table 4 indicates that the maximum dis- 
tance probed by detectors of first generation (VIRGO, 
LIGO I) is about 6.4 Mpc, well beyond M31, whereas 
the second generation (LIGO II) will probably see phase 
transition events at distances two times longer, but cer- 
tainly not yet attaining the Virgo cluster. The small 
probed volume and the rapid rotation required for this 
mechanism be efficient imply in a low event rate, im- 
posing severe limitations on the detectability of such a 
signal. 
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